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Random free fermions: An analytical example of eigenstate thermalization. 
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Having analytical instances of the Eigenstate Thermalization Hypothesis (ETH) is of obvious 
interest, both for fundamental and applied reasons. This is generically a hard task, due to the 
belief that non-linear interactions are basic ingredients of the thermalization mechanism. In this 
article we proof that random gaussian free fermions satisfy ETH in the multiparticle sector, by 
analytically computing the correlations and entanglement entropies of the theory. With the explicit 
construction at hand, we finally comment on the differences between fully random Hamiltonians and 
random Gaussian systems, providing a physically motivated notion of randomness of the microscopic 
quantum state. 
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I. INTRODUCTION 

The problem of quantum thermalization can be stated 
as: 

• Given microscopic unitarity, how do Gibbs ensem¬ 
bles emerge? 

If a many body quantum system is set in an initial pure 
state |i/’(0)), the evolved state = U{t)\ip{0)) is pure 
and time dependent, so it can never become a time inde¬ 
pendent mixed density matrix, such as Gibbs ensembles: 

( 1 ) 

where is any Gibbs distribution. The dynami¬ 

cal emergence of Gibbs ensembles from unitary dynamics 
will be termed the problem of quantum thermalization. 
This problem is almost as old as quantum mechanics it¬ 
self. For a selfcontained recent review, with an excellent 
account of references and historical rigor, see [T|. 

Although exact thermality cannot be attained within 
unitary evolution, we might still expect approximate 
thermality for the actual measurements done in exper¬ 
iments, the correlation functions of the theory. Given an 
observable O of the theory, it is written mathematically 
as: 

(V’(i) \0\^{t)) = Tr(p°'‘’'"’^ O) ± error . (2) 

stating that the correct expectation value, as measured 
by the evolving quantum state is equal to a Gibbs en¬ 
semble average, up to some error. This error has to be 
negligible in the thermodynamic limit for the previous re¬ 
lation to be non trivial, and for the thermal expectation 
value to be a good approximation of the correct one. 

To understand such behavior it was proposed in Q] 
that the typical energy eigenstate \Ea) of the quantum 
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system satisfy themselves the previous property: 

{Ea\0\Ea)= Tr(pGibbs ± ^ (3) 

stating that the expectation value in the energy basis is 
well approximated by the thermal correlator. The previ- 
ous phenomenom was coined Eigenstate Themalization 
in [ij, but was earlier noticed in Q, in the context of 
quantum spin systems. At any rate, as argued in [i|, 
both previous equations can be seen as more precise for¬ 
mulations of the original question (|T|, and as such they 
are hypothesis about the nature of quantum systems. In¬ 
deed, the second equation (jS]) is commonly known as the 
Eigenstate Thermalzation Hypothesis (ETH). 

Relations m and m suggest various obvious ques¬ 
tions: 

• What type of systems display such behavior? Over 
which range of initial states and eigenstates do they 
display it? 

• What is pGibbsv ^]^g effective temper¬ 

ature T = 1//3 from the pure state scenario? 

• What is the ‘error’? 

Although we have a great deal of intuition about these 
questions through the connection between ETH and ran¬ 
dom matrices, to the author knowledge, see [l|, no time 
independent Hamiltonian containing only few body inter¬ 
actions has been proven to satisfy relations @ and ©■ 
The difficulty to prove such behavior is related to the be¬ 
lief that non-linear interactions are fundamental to the 
quantum thermalization mechanism. 

The objective of this article is to study a family of 
quantum systems displaying ETH, satisfying relation ([3]). 
Gontrary to common belief, these systems are ‘gaussian’, 
containing only two-body interactions, see (|1]). We fur¬ 
ther compute the errors, both in correlations and entan¬ 
glement entropies. These errors differ from the random 
matrix approximation, providing a physically motivated 
measure of ‘randomness’ of the quantum state. 
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We want to remark that although these systems might 
seem unnatural from several points of view, they indeed 
serve as microscopic toy models of black physics. The 
underlying reason is that they are examples of systems 
with no locality structure whatsoever. The Hamiltonian 
connects every oscillator with every other democratically, 
a characteristic feature of black hole physics 

From a different perspective, since these systems dis¬ 
play ETH, entanglement entropy is extensive across all 
bipartitions, even for ‘vacuum states’ in the sectors with 
a large number of particles. This could be of interest for 
the type of questions raised in 0 . 

II. RANDOM FREE FERMIONS, 
ENTANGLEMENT ENTROPY AND ETH 

The family of Hamiltonians we wish to study is the 
following: 

N N 

H = a'^cl c, + r] cl V,j cj , (4) 

i—1 i,j—l 

where a and ry are parameters with energy dimensions, cj 
and Ci are creation and annihilation operators of spinless 
free fermions, with usual anticommutator relations, and 
the couplings Vij are independent random gaussian real 
numbers with zero mean and unit variance. The matrix 
{r]V)ij = r]Vij is therefore a random matrix taken from 
the GOE ensemble with deviation ay = V (see [l^ for a 
beautiful and modern treatment of random matrices). 

The ‘free’ nature of the model allows an exact so¬ 
lution via diagonalization of the matrix V. If ^/>“, for 
a = 1, • • • ,N, are the eigenvectors of V with eigenvalues 

Ca- 

N 

= ( 5 ) 

then the Hamiltonian can be written as: 

N N 

H = '^{a + ea)dlda = '^Eadlda , (6) 

a—1 a=l 

where d^ and da are new creation and annihilation oper¬ 
ators defined by: 

N 

da = Y^ Ci . (7) 

i=l 

All eigenstates are constructed by choosing a set A of 
particles, and associate to it the following eigenstate: 

n^aio), (8) 

aeA 

where | 0 ) is the state annihilated by all da, and Np is the 
number of particles in the state, a number which will play 


a key role below. Notice that there are ) independent 
states for a given Np. 

The objective of this article is to study the structure of 
these eigenstates, associated with 0 . We will focus on 
the number operator, the two point correlation functions 
and entanglement entropies. Since the correlations and 
entanglement structure of the eigenstates are symmetri¬ 
cal under Np -A N — Np , a manifestation of particle-hole 
symmetry in this model, we will only focus on Np < N/2. 

To compute the correlations and entanglement en¬ 
tropies, we will make use of the theory of random matri¬ 
ces, which deals with the statistical properties of eigenval¬ 
ues and eigenvectors of matrices such as rjV (see [T^ , [3 
and poj i for an extense treatment of random matrices). 
In relation to the eigenvectors, the main assertion is that 
the orthogonal matrix of eigenvectors ,■■■ , ip^) is dis¬ 
tributed according to the Haar measure on the orthog¬ 
onal group 0{N). For our purposes, this means that 
the eigenvectors have independent and random gaussian 
entries, up to normalization. Matemathically: 

K] = 0 bPt ^j] = Sab , (9) 

where [p] denotes the average of the random variable p 
over the matrix ensemble. 

For the eigenvalues we will only need Wigner’s semi¬ 
circle law, accounting for the probability of having an 
eigenvalue equal to A: 

^ ’ (10) 

where = A:Nrf, and where we remind that this law 
concerns the eigenvalues of 77 H, a matrix with deviation 
equal to 77, see the Hamiltonian (|4]). To ensure a zero en¬ 
ergy vacuum eigenstate of H we assume R = 2\fNr] < a. 
Given CQl), the first two moments for the eigenenergies 
Ea are: 

[Y,Ea]=Npa [(^K)2]-[^Kp=iV,iV77^ 
a^A a^A a^A 

( 11 ) 

With the previous statistical information about eigenval¬ 
ues and eigenvectors we can now compute the correlation 
functions for the state ([5]) with Np particles: 

N 

= (4/^-1 clc, 1 4 /^-)= ^ (4/^-14*14-^-) = 

a, 6 — 1 

= ( 12 ) 

a£A 

where the last sum just runs over the subset A of par¬ 
ticles chosen. The previous correlator is itself a random 
variable, a functional of the random eigenvectors 7/1“. By 
using we obtain: 

[C?,] = ^ S,, [Cf.Ci] - [Cfj] [Cf,] = ^ 5,u Sp . 

( 13 ) 
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We thus can see the correlation matrix as the sum of 
a ‘thermal’ part and an error, the error being a ran¬ 
dom matrix taken from the GOE ensemble with deviation 
Oq = ^. Notice that the ‘thermal’ part just depends on 
the ‘macroscopic’ parameter Np, and cannot distinguish 
between the ) independent states in the correspond¬ 
ing Np sector. Also notice that the ‘thermal’ part is a 
good approximation in the window Np ^1. As might 
have been expected, and anticipating results, Np = N/2 
will correspond to the high temperature sector. 

To proof ETH for this system, relation ([3]) described in 
the introduction, we need to compute the corresponding 
correlation matrix in the thermal ensemble. This is given 
by: 

= |Tr(e-^^clc,) 

N 

a—1 

where = l/(e^'®“ -|- 1) is the average number oper¬ 
ator for a fermionic oscillator at temperature T = 1//3, 
and Z = Tr(e“^^) is the usual partition function. The 
thermal correlation matrix is again a random variable, 
due to the randomness of H. Since there is no correla¬ 
tion between eigenstates and eigenvalues, the mean and 
variance are given by: 

[C,§] = [np] % - [C,§] [Cfj = <5,, . 

(15) 

where: 

R R 

[np] = I P{X)4d\ [inpf]= I P{X)in^pfdX. 

-R -R 

( 16 ) 

We did not succeed in analytically computing the aver¬ 
ages for all /3. In the large temperature limit. Pa —>■ 0, 
they read: 

whereas in the low temperature limit, PR —>■ oo, we have 
slightly more complicated expressions: 


we conclude that the leading diagonal approximation of 
the thermal correlation matrix is valid at all tempera¬ 
tures. 

Now we are ready to compare both correlation ma¬ 
trices, the exact eigenstate correlations m versus the 
thermal (USD ones, allowing us to arrive to the following 
conclusions: 

• For Np ^ 1, a pure eigenstate with Np particles 
can be effectively approximated by a Gibbs distri¬ 
bution at temperature P satisfying ^ = [np], the 
difference being subleading in the thermodynamic 
limit. Random free fermions then constitute an ex¬ 
plicit analytical example of ETH (jSD, an example in 
which we know the effective temperature and the 
error size in terms of the microscopic parameters of 
the theory. 

• For Np ^ 0{1), approximating pure states by ther¬ 
mal ensembles is not a valid approximation. 

These features can be made clearer by studying entan¬ 
glement entropy. For gaussian systems, as shown in [^ . 
one can compute the entanglement entropy of a given 
subsystem A directly from the correlation matrix. More 
concretely, given a subsystem A with m degrees of free¬ 
dom, knowledge of = ('I'|cJcjjd>), where i,j G A, 
allows for the computation of the entanglement entropy. 
It is given by: 

m 

5A = -^(AlogA, + (I-A)log(I-A)) , (20) 

i=l 

where Aj, with i = 1, ■■ ■ ,m are the m eigenvalues of the 
matrix in the given subsystem A. The proof relies 
only on the fact that the correlation matrix and the re¬ 
duced density matrix share the same set of eigenvectors. 

Although formula (|20l) is fairly simple, one still need 
to compute the eigenvalues of , and this is not always 
possible analytically, even for one-dimensional systems, 
see [2l|. Indeed we were not able to compute the entan¬ 
glement entropy for all Np. We will compute it in the 
two standard limiting cases: the thermal regime, speci¬ 
fied hy Np 1, and for Np = 1, corresponding to the 
non-thermal phase. 

Let us begin with the simpler case of having just one 
particle. The wave function is then given by: 


= ^ Tr(e 4 4) = 

a, 6 — 1 

(14) 


[np] = 


2e-“^Ji(/3i?) 

PR 


0-/3 (a-R) 


nPR PR 


and 


[inpf] = 


g-2/3a g-2/3(a-iJ) 


PR 


PRV^^ ’ 


(18) 


(19) 


where Ii{x) is the modified Bessel function of the first 
kind. Since to ensure a zero energy vacuum eigenstate of 
p[ we assumed R = 2y/Nf] < a, from (|T5D . (|TOD and (fTOl) 


N N 

= 4 |o) = E 4 |o) = E V’? 4 > (21) 

where |i) = c| |0). This state is considered in in 
relation to the Many-Body-Localized phase transition. 
Because the state is fully supported in the single parti¬ 
cle sector, the entanglement of any subsystem A can be 
expressed as: 

Sa = -PA logPA - (1 - pa) log(l - pa) , (22) 
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where pa = X] see [l^ for an explicit derivation of 

ieA 

the previous formula. The average and variance of pA are 
[pa] = m/N and 0 -^^ = m/N'^. Because the variance of 
the probabilities pA is small in comparison to the mean 
value, we Taylor expand to compute the average entan¬ 
glement entropy, following Q for the case of random QFT 
states. For the case at hand, the entanglement entropy 
of a subsystem A oi m < N/2 degrees of freedom in the 
state m is given by: 


[m] N ^ 2(lV-m)’ 

(23) 

where we remark that on average does not depend on 
the particle type chosen, labeled by a, and that the last 
term in the previous equation is alway subleading in 
the thermodynamic limit. The previous formula for en¬ 
tanglement entropy is a monotonic growing function for 
1 < TO < iV/2, with S'?' ~ and ~ log 2. This 
is a very small amount of entanglement, and certainly 
not extensive with the number of sites to. Therefore it 
cannot be faithfully represented as a thermal entropy, 
signalling that the small Np sector does not satisfy ETH. 

For the high Np sector we follow a novel approximate 
method developed in Q , Appendix A. The key aspect to 
observe is that for Np ^ 1 the correlator matrix m is 
a diagonal matrix plus a random matrix with paramet¬ 
rically smaller entries. Writing = C'^ + 6C'^, to first 
order the eigenvalues of are given by Xi ~ Xi + SXi, 
where Xi = Np/N and 6Xi are the eigenvalues of a ran¬ 
dom matrix of size to with deviation , see (ITSl) . 

Using Wigner’s semicircle law (nni, for a matrix of size to 
with such variance, these eigenvalues satisfy [i5Ai] = 0 and 
[((5Ai)^] = TO^ . Therefore, for m < Np, plugging these 
eigenvalues into (l20l) . Taylor expanding in 6Xi and finally 
taking the average, we obtain the following formula for 
the average entanglement entropy of a subsystem of size 
TO < Np in the multiparticle state 




m 


2(iV - Np) 


where we have defined: 


^ ^ _ (1 _ ^) log(l _ 

1 N ^ N ^ N ^ N ’ 


(24) 


(25) 


corresponding to the thermal entropy per degree of free¬ 
dom for a state with Np ^ 1 particles. Notice that to 
change from Np to the effective temperature T = 1//3 we 
just need to use ^ = [np]. 

Given the explicit expressions (IMl) and (EH) we con¬ 
clude: 


• Entanglement entropy is well approximated by the 
thermal entropy for subsystems with sizes smaller 
than the number of particles in the given eigen¬ 
state. For bigger subsystems we cannot state any¬ 
thing with certainty, but it is tempting to specu¬ 
late a slower growth of entanglement entropy for 


^ ^p- For Np N/2 the thermal approxima¬ 
tion is valid for every subsystem. 

• The random nature of the state seems to increase 
with the number of particles. 


III. ERROR SCALINGS AS RANDOMNESS 
MEASURES 

Until now we have focused in the structural properties 
of the eigenstates of the Hamiltonian (U), in particular 
on their thermal nature. It is now fruitful to check the 
differences between traditional approaches to quantum 
chaos and ETH based on fully random Hamiltonians, and 
the exact solution we found. This will provide a physi¬ 
cally motivated quantification of the ‘randomness’ of the 
quantum state. 

Although this section is self-contained, for extensive re¬ 
views about the traditional approach to ETH and quan¬ 
tum thermalization based on ‘typicality’ arguments and 
random matrices in the contexts of quantum informa¬ 
tion theory, condensed matter and black hole physics see 
[liQ- To proceed we first write the Hamiltonian (U) 
as: 


H = Hfree + V Hint ■ (26) 

In most cases in which Hint contains non-linear inter¬ 
actions we cannot diagonalize H exactly. The tradi¬ 
tional approximation assumes that the Hamiltonian can 
be taken from one of the random matrix ensembles, see 
. The approximation implies the famous chaotic spec¬ 
tra, and more concretely Wigner’s semicircle law m- 
More interestingly for the concerns of this article, it also 
implies ETH as we show below, see 

Let us begin with the spectrum. Since the dimension 
of the subspace with Np particles is (^), and assuming 
the entries of Hint to be random gaussian variables with 
zero mean and unit variance, the eigenenergies satisfy: 


Er‘^°^ = Npa±2^(^^^fj, (27) 

where a = I,-- - , {^)- The first two moments of the 
random approximation can be computed from Wigner’s 
semicircle law: 

j^random] ^ ]S[p a [(£;™dom^2] _ j^random]2 ^ ^ -2 

^28) 

On the eigenvectors side, the assumption of a random 
Hamiltonian implies that the eigenvectors are random 
vectors in the corresponding subspace: 

l^random) ^ g ^ (29) 
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where |i) is a basis in the Np subspace and: 

= 0 [C ilJj] = 6ab ■ (30) 

\Np) 

See [1] for a recent detailed treatment of these type of 
vectors. Using dSHD, the statistical properties of Clj = 

^^randomj j^random^ 


[cy = 


Hc'.if] - [cijf = 


/Ar-2N 

V7Vp-l/ 

\nJ 


OiTwr)- 

\nJ 


(31) 

Comparing formulas (1^51) and (EU, corresponding to the 
random Hamiltonian approximation, with the exact so¬ 
lutions (HU and m, and with the thermal result (HU, 
we conclude that: 


• Appropriately fixing fj in terms of ry, Np and TV, the 
mean and variance of the eigenvalues is the same 
for the exact Hamiltonian (|3]) than for its random 
approximation. 

• The average correlation matrix [Cij] coincides for 
all cases, the differences lying on variances. These 
variances furnish good quantifiers of randomness 
in the quantum state. They subtly distinguish be¬ 
tween macroscopically equal phases, such as eigen¬ 
states of fully random hamiltonians and eigenstates 
of random free fermions. The scaling properties of 
the errors in typical eigenstates seem a fruitful field 
to explore in the context of quantum thermaliza- 
tion. 

• The random Hamiltonian approximation seems 
valid for all subsystem sizes m (notice that in the 
one particle sector is exact for all subsystems), 
while the Gibbs distribution seems to hold only for 
m < Np. 


IV. CONCLUSIONS 

In this article we studied aspects of the Hamilto¬ 
nian O, such as the spectral properties (HU), correlation 
matrix HU and entanglement entropies, (HU and (HU- 
They can be expanded in 1/iV, where N is the number of 
spinless fermions of the model. The leading term in this 
expansion is always the thermal result, given by HU, 
the effective temperature being found in terms of the 
macroscopic parameter Np characterizing the sector. The 


family of Hamiltonians 0 thus satisfy 0, furnishing 
explicit examples of the Eigenstate Thermalization Hy¬ 
pothesis (ETH). The conclusion is remarkable, since it 
implies that ETH is typical within the space of gaussian 
Hamiltonians. Since it is also typical in the full space of 
Hamiltonians [I| , as proved by relation m, it is tempt¬ 
ing to conclude that it is typical for Hamiltonians with 
random 2,3 •••TV body interactions, such as the model 
presented in fl^ . 

We found that the entanglement properties of big sub¬ 
systems in these eigenstates are different from the ther¬ 
mal result. For a sector with Np ^ 1 particles, we were 
able to prove thermality of entanglement entropies un¬ 
til subsystems of size of 0{Np). For bigger subsytems we 
speculated with a slower growth of entanglement entropy, 
but otherwise further study is needed to unravelled its 
nature. These results might have impact in black hole 
physics, in which the present picture is that given by 
Page in [2^, a picture in which thermality holds for ev¬ 
ery subsystem, and in which deviations from thermality 
are assumed to be those given by the random Hamilto¬ 
nian approximation. 

The last section was devoted to study the differences 
between the approach to quantum chaos and ETH based 
on random matrices, see the reviews Bin , and the 
exact solution of the Hamiltonian (H]). These differences 
lie in the deviations from the thermal result. The ‘errors’ 
in equations dU and ([U are good quantifiers of the ‘ran¬ 
domness’ of the quantum state. Their scaling properties 
seem an exciting new route to explore in the context of 
quantum thermalization. 
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